Softly Constrained Films 
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The shape of materials is often subject to a number of geometric constraints that limit the size 
of the system or fix the structure of its boundary. In soft and biological materials, however, these 
constraints are not always hard, but are due to other physical mechanisms that affect the overall force 
balance. A capillary film spanning a flexible piece of wire or a cell anchored to a compliant substrate 
by mean of adhesive contacts are examples of these softly constrained systems in the macroscopic 
and microscopic world. In this article I review some of the important mathematical and physical 
developments that contributed to our understanding of shape formation in softly constrained films 
and their recent application to the mechanics of adherent cells. 



I. INTRODUCTION 

A typical problem in the mechanics of soft materi- 
als consists in finding the optimal shape of a system in 
the liquid phase subject to a number of geometric con- 
straints. The problem occurs on the macroscopic scale, 
when one would like to determine the shape of a droplet 
wetting a surface [1] or that of a soap film spanning a rigid 
wire frame [2]. In the microscopic world, the mechanics 
of membranes, vesicles or micro-droplets revolves around 
finding the shape the minimizes the energy of the sys- 
tem given its volume [3, 4]. That on the volume (or the 
area in two dimensions) is a hard constraint, that results 
from the incompressibility of matter in the liquid phase. 
Similarly, imposing the shape of the boundary spanned 
by a two-dimensional interface, leads to vast class of con- 
strained shape optimization problems in capillarity [5]. 
It is then natural to ask how the problem would changed 
if the volume, area or boundary constraint was relaxed. 

Soap films bounded by elastic rods provide a table top 
realization of these softly constrained systems. A classic 
demonstration, when familiarizing with the concepts of 
surface tension, consists of dipping a rigid wire frame into 
soapy water and let the soap film form a surface of zero 
mean curvature (i.e. a minimal surface) [2]. How would 
the shape of the film change if the bounding frame was 
not rigid, but rather made of an inextensible but flexible 
elastic rod? And what if the boundary was made of a vis- 
coelastic layer whose length itself was softly constrained 
by other physical mechanisms such as adhesion? 

The first example of this class of problems dates back 
to the end of the 19th century in the work of Maurice 
Levy [6], who investigated the stability of an infinitely 
long cylindrical pipe under pressure. Because the sys- 
tem is invariant for translation along the longitudinal 
direction, all the forces act on the cross-sectional plane 
of the pipe and the system is formally equivalent to a 
two-dimensional capillary film enclosed by a flexible elas- 
tic rod. This problem drew the attention of many re- 
searchers for more than a century and was finally solved 



a few years ago in a completely analytical manner [7, 8] . 
Although the physics involved is very simple and ulti- 
mately reduces to the competition between two antag- 
onist forces that try to bend the boundary in opposite 
directions, the problem provides a tremendous example 
of polymorphism and multi-stability. By simply increas- 
ing a single physical parameter, the surface tension, one 
can access an increasingly large number of stable config- 
urations and in the limit of infinite surface tension the 
system admits a countably infinite number of equilibria. 

The Euler-Plateau problem, is a generalization of 
Levy's original problem obtained if the softly constrained 
film is allowed to bend in three dimensions. The prob- 
lem was proposed in Ref. [9] as a link between two of 
the oldest problems in geometrical physics and the cal- 
culus of variations, bringing together Euler's Elastica and 
Plateau's problem. The union of these leads to the ques- 
tion of the shape of minimal surfaces bounded by elastic 
lines, a new class of questions at the nexus of geometry 
and physics. Among biological systems the paradigm of 
softly constrained film was applied in tissue morphogen- 
esis to account for the looping patterns in the developing 
vertebrate gut tube [10] and in various two-dimensional 
models of adhering cells [11-13]. 

This short review is organized as follows. In Sec. II 
we will refresh some basic concepts of differential geom- 
etry of curves, mostly to establish notation. In Sec. Ill, 
we will consider the classic problem of a two-dimensional 
film bounded by an elastic rod, we will derive the ana- 
lytical solution and discuss in detail the various regimes 
obtained by varying the surface tension of the film. The 
Euler-Plateau problem for softly constrained soap films 
in three dimensions will be discussed in Sec. IV. In Sec. 
V we will focus on the Contractile Film Model for adher- 
ent cells [13] and see how the notion of softly constrained 
films finds a natural and important application to cell 
mechanics. Conclusions will be drawn in Sec. VI. 



II. MATHEMATICAL PRELIMINARIES AND 
NOTATION 
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In all the upcoming sections, we will make a large use 
of the Frenet-Serret apparatus for curves in two and three 



dimensions, whether immersed in Euclidean space or ly- 
ing on a surface. In this section we will briefly review 
some fundamental concepts, mostly to establish nota- 
tion. Let r(s) be the position vector of a curve in R 3 
parametrized by its arc-length s. The Frenet frame of 
the tangent t = dr/ds, normal n and binormal vector b 
is described by the Frenet-Serret formulas: 



i! = nn . 



n 



-Kt + rb , 



(la) 
(lb) 



6 / = -r6, (lc) 

as well as the following cross-product relations: 



t = n x b 



n 



bxt 



t x n . 



(2) 



Here k and r represents respectively the curvature and 
the torsion of the curve and ( • )' = d(-)/ds. Eqs. (1) 
describe how the orthonormal frame {£, n, 6} rotates in 
space as we move along the curve (Fig. 1, left). Once n 
and r are assigned, in the form of two different iable func- 
tions, the fundamental theorem of space curves guaran- 
tees that the corresponding curve is uniquely determined, 
up to a rigid motion [14]. Note that in three dimensions, 
the equation t' = nn does not define the sign of the 
curvature. That is conventionally chosen to be always 
positive, by incorporating the sign of t' in the definition 
of n: i.e. n = t' /\t'\. With this choice, the normal vector 
n is oriented toward the center of curvature of the curve 
(i.e. the center of the osculating circle). It is useful to re- 
call that the Frenet frame is well defined only where the 
curvature n does not vanish. Where n = 0, the normal 
vector n = t' j k is undefined and so are b and r. 

If the curve lies on a plane (Fig. 1, right), its tan- 
gent vector can be conveniently parameterized through a 
single scalar function 0(s) representing the turning angle 
of the tangent vector t. In a standard Cartesian frame: 
t = (cos#, sin#), then, using Eq. (1), one finds n = 9' 
and r = 0. In this case, the curvature n identifies the 
curve uniquely, up to rigid motion: 

r(s) = r (0) + / ds' [cos 0(s') x + sin 0(s') y] (3) 
Jo 

with 0(s) = 0(0) + J ds' k,(s). Expressing the curvature 
of a plane curve as the derivative of the turning angle 
0, allows one to unambiguously identify the sign of n: 
k > (k < 0) implies that the turning angle increases 
(decreases) as we move along the curve, while reversing 
the orientation of the curve changes the sign of n. If a 
closed plane curve is oriented counterclockwise and the 
turning angle is measured, as usual, from the x— axis of 
a Cartesian frame, k, > (k, < 0) will then corresponds 
to points where the curve is convex (concave). This con- 
vection will be used throughout the paper. 

Closed curves, have a number of interesting global 
properties that will serve as important calculation tools 
in the later sections. The four vertex theorem [15], is 




FIG. 1: The Frenet frame in three (left) and two (right) 
dimensions. In three dimensions the curvature k is conven- 
tionally taken with the positive sign and the normal vector n 
is directed toward the center of curvature. If a closed plane 
curve is oriented counterclockwise, on the other hand, n is 
directed toward the interior of the curve and k > (k < 0) 
corresponds to points where the curve is convex (concave). 



one of the earliest results in global differential geometry; 
it states that the curvature of a simple, smooth, closed 
curve on the plane has a least four vertices: i.e. four 
extrema where k' = (specifically two maxima and two 
minima). Another fundamental property of closed plane 
curves is expressed by the theorem of turning tangent 
[14, 16]. This states that: 



ds k = 27rm . 



(4) 



The integer m is called the rotation index of the curve 
and measures how many times the curve turns with re- 
spect to a fixed direction [16]. Simple closed curves have 
thus m = 1 (Fig. 2, left), while a curve that loops twice 
around its center (thus self-intersects once before closing) 
has m = 2 (Fig. 2, center). If a simple closed curve has 
kinks (singular points where the tangent vector switches 
discontinuously between two orientations), these will re- 
sult in the total curvature as it follows: 



ds n 



£< 



2tt 



(5) 



where fa is the external angle at each kink and the sum- 
mation runs over all the kinks (Fig. 2, right). In the case 
of a convex polygon, for instance, k = and (5) asserts 
that the sum of the external angles is equal to 2ir. In the 
case of a closed space curve, a theorem by Fenchel [17] 
states that: 



ds k, > 2tt 



(6) 



where the equality holds only if the curve lies on the 
plane. 

A space curve that lies on a curved surface carries an 
important information about the local geometry its non- 
Euclidean habitat. This can be highlighted by introduc- 
ing the Darboux frame consisting of the curve tangent 




FIG. 2: Three closed plane curves and their total turning 
angle AO = § ds k. (Left) In a simple closed curve, like a 
circle, the tangent vector rotates by 2ty in one loop around 
the curve. (Center) In the case of a limacon, the tangent 
makes two full rotations, thus AO = 4tv. (Right) A deltoid 
curve has three cusps where the tangent vector jumps by an 
angle <j) — tc. Eq. (5) then requires § ds n + 3tt = 2tt, thus 

AO = -7T. 



vector £, the surface normal vector N and the so called 
tangent-normal vector (i.e. tangent to the surface and 
normal to the curve) T = N x t [18]. The Darboux and 
Frenet frame can be transformed into each other through 
a rotation by an angle = Z(6, N) about t that brings 
b onto N and n onto T. In matrix form: 



t \ 1 

T = cos 
N \ -sin 




(7) 



Then, taking a derivative with respect to the arclength s 
and using the Frenet-Serret formulas (1) one obtains the 
equations: 

t' 



T' = 

N' 



Kg J- \ K" 


JV, 


-- —Kgt + 


r g N, 


= —K n t — 


r g T. 



(8b) 
(8c) 



that describe how the orthonormal frame {£, T, N} ro- 
tates as we move along the curve. n g = ftcos</> is the 
geodesic curvature, n n = — ftshi(/> is the normal curva- 
ture, while r g = r + <j>' is the geodeisc torsion. The 
geometrical meaning of k q and K n becomes evident by 
writing Eq. (8a) in the form: 



„T + K n N . 



(9) 



K n and Kg are the projections of the curvature vector ku 
along the surface normal and tangent direction respec- 
tively. Thus n n represents the contribution to the curva- 
ture due to the fact that the underlying surface is itself 
curved, while n g is the intrinsic curvature of the curve. 
Geodesies, being the most direct path between two points 



on a surface, have K q — 0, while curves with n n 







are called asymptotic lines. An alternative parametriza- 
tion of the Darboux frame (which we will adopt in Sec. 
IV), makes use of the angle $ = + tt/2 between the 
curve normal n and the surface normal N. In this case: 
Kg = ft sin $ and n n = n cos r d. 



III. PLANAR FILMS BOUNDED BY ELASTIC 
RODS 

In Sec. I we introduced the problem of a soap-film 
bounded by an elastic rod (Fig. 3) as the prototype 
problem for softly constrained films. Here we review the 
simpler (and yet incredibly rich) version of the problem, 
in which the system is constraint to lie on a plane and its 
geometry is entirely determined by that of the bound- 
ary curve. The problem has a longstanding history as 
it maps one-to-one to the important problem of finding 
the cross sectional shape of an infinitely long cylindri- 
cal pipe subject to a pressure difference. Since the pipe 
is flat along its longitudinal direction, its bending en- 
ergy is completely determined by the curvature of the 
boundary and the system is physically equivalent to a 
two-dimensional rod (representing the boundary of the 
pipe) subject to a constant normal force (i.e. the pres- 
sure difference). 

Levy [6] was the first one to formulate the problem at 
the end of the 19th century. He investigated the stabil- 
ity of the circular shape and reduced the calculation of 
the non-circular solutions after buckling to an algebraic 
problem involving elliptic integrals. Later, Halphen [20] 
and Greenhill [21] derived an exact solution in terms of 
Weierstrass elliptic functions and unveiled the existence 
of a sequence of buckling transitions to shapes of increas- 
ing order in rotational symmetry for increasing surface 
tension. This solution was however considered too com- 
plicated and unintuitive by many researchers who con- 
tinued looking for alternative exact and approximated 
solutions. 

Some fifty years after the original work by Levy, Car- 
rier resumed the problem, found an expression for the 
curvature of the buckled ring in terms of Jacobi elliptic 
functions and worked out an approximate solution near 
the onset of buckling. This analysis as been recently 
reconsidered by Adams [22] who further extended Car- 
rier's work. Tadjbakhsh and Odeh [23], provided a rigor- 
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FIG. 3: A rod enclosing a film of a commercial dishwashing 
liquid. The Young modulus of the rod is E — 3.55 ±0.1 GPa, 
the radius is p — 0.25 mm, and the surface tension of the 
film is a = 27.3 mN/m. The dashed curve is obtained from 
the theoretical solution of Eq. (34) with no fitting parameter 
and is described in detail in Sec. IIIC. Reproduced from Ref. 
[19]. 



ous study of the boundary- value problem describing the 
shape of the ring and the associated variational problem. 
In addition, they performed the first numerical analysis 
of the problem and visualized a variety of post-buckling 
solutions, including complicated self- intersecting curves. 

Flaherty et al. [24] lifted the problem to a whole new 
level by proposing a novel scenario for the evolution of 
the equilibrium shape for increasing surface tension that 
includes self-contact. According to this picture, upon in- 
creasing the surface tension above a threshold, the ring 
collapses into itself. The collapse initially occurs with the 
appearance of an isolated point of contact, that for larger 
surface tension expands into a segment, giving rise to a 
shape reminiscent of a two-headed tennis racket. Inter- 
esting physical and mathematical aspects of the problem 
have been further investigated by Arreaga et al. [25], 
Capovilla et al. [26] and Guven [27]. 

A completely analytical solution was finally given in a 
remarkable series of papers by Vassilev, Djondjorov and 
Mladenov for shapes before [7] and after collapse [8] . The 
latter case has been recently revisited by Mora et al. [19] 
who additionally provided and insightful comparison be- 
tween analytics and experiments (Fig. 3). In the fol- 
lowing section, we will review some fundamental aspects 
of the solution and highlight some qualitative features, 
while remanding to the original works for a fully detailed 
treatment. 



A. Shape equation 

Let us consider the following functional representing 
the total mechanical energy of the system: 



E 



a dA+ (b i 

JM JdM 



ds (an 2 + /3) 



(10) 



Here M is a mapping of a two-dimensional disk represent- 
ing a film of surface tension <r, constrained to lie of the 
plane and bounded by a flexible beam of circular cross- 
section with bending rigidity 2a. The boundary dM is 
thus treated as an Elastica immersed in two-dimensional 
Euclidean space 1R 2 , and /3 is a Lagrange multiplier, anal- 
ogous to line tension, that enforces inextensibility of the 
boundary: 



L 



ds , 



(ii) 



dM 



with L the perimeter. For the two-dimensional problem 
discussed here, we used the convention introduced in Sec. 
I, thus the tangent vector t is oriented counterclockwise, 
n is directed toward the interior of the loop, so that the 
curve is convex where n > and concave where k, < 
0. As a mechanical energy, Eq. (10) is amenable to 
various interpretations. Thinking of a as a hydrostatic 
pressure, rather than a surface tension, and to k as the 
only nonzero principal curvature of a cylinder obtained 
by extruding dM along the z— axis, Eq. (10) is equivalent 



to the energy of infinitely long cylindrical pipe subject to 
a uniform hydrostatic pressure. Analogously, replacing 
the constraint on the length (11) with a constraint on the 
area, with a playing the role of a Lagrange multiplier, and 
interpreting /3 as a tunable line tension, Eq. (10) becomes 
the energy of a two-dimensional vesicle consistently with 
the classic framework of Canham [28] and Helfrich [29]. 
To derive the conditions for the extremum of the func- 
tional (10), we will impose a small virtual deformation of 
the boundary, r — >■ r + £n, and use this to calculate the 
corresponding change in energy and the necessary condi- 
tions for the resulting shape to be a local minimum. A 
linear variation of the boundary energy can be calculated 
with standard methods in the following form [30]: 



5 I ds {an 2, 

JdM 



'P) 



ds [a(2n" + n 3 ) - /3«]f (12) 



dM 



The contribution due to surface tension, can also be read- 
ily found by expressing the enclosed area as a contour in- 
tegral: A = —(1/2) § dsr • n. Then, taking into account 
that S(ds) = —£K,ds and S(n) = — £'£, the area variation 
can be readily found to be: 



SA = 



ds£ 



(13) 



dM 



consistently with intuition. Combining Eqs. (12) and 
(13) and setting the integrand to zero, we finally obtained 
the Euler-Lagrange equation: 



ql(2k" + k 3 ) - /3k - a = 



(14) 



Physically Eq. (14) represents the balance between the 
normal component of the stress resultant F acting a 
cross-section of the rod and the body force K = —an 
due to surface tension. Namely: 



F' = an. 



(15) 



In two dimensions, the stress resultant F can be found 
using the principle of virtual works [9, 27]: 



F = (an 2 - 0)t + 2<Wn 



(16) 



Differentiating Eq. (16) and using Eq. (15) readily gives 
Eq. (14). 

Eq. (14) belongs in the class of equations describing 
conservative systems with one degree of freedom [7] . In- 
tegrating Eq. (14) with respect to ft, gives: 



E = a(K') 2 + U , 
where U is the fourth-order polynomial: 



(17) 



U = \ol^ - \j3n 2 - gk , (18) 

and £ an integration constant. Thus, interpreting the 
quantity s as a time and n as a position, Eq. (14) corre- 
sponds to the equation of motion of a particle of mass 2a 



and kinetic energy a(n') 2 , moving in a potential U and 
whose total energy £ is conserved. 

Important geometric information about the shape of 
this two-dimensional softly constraint film, can be found 
directly from Eqs. (15) even before integrating Eq. (14). 
This was first established in [25] and later expanded in 
[26, 27]. It is not hard to verify that the normal vector n 
can be expressed as an arc-length derivative of a vector 
field: 



n = -j-[(r • t)n — (r • n)t] . 



(19) 



Thus, the force balance relation (15) can be recast as an 
identity for the plane vector: 

[an 2 - j3 + a(r ■ n)]t + [2<W - a(r ■ t)]n = C , (20) 

where C is some constant vector along the loop. This 
constant is unimportant as one can always choose the 
reference frame in such a way that r • t — and r • n — 
— {an 2 — /3)/<j where n' = 0, which implies C = 0. For 
a 7^ 0, this allows one to parametrize the curve in the 
Frenet frame it follows: 



F x z 2an' 



an^ 



p 



a 



a 



(21) 



where the first equality is obtained from Eq. (2) by set- 
ting z = b. Because of the symmetry of the equilibrium 
shapes, § ds F = 0, which implies that the origin of the 
reference frame associated with Eq. (21) coincides with 
the center of mass of the curve [26]. Using this and (17), 
one can further derive a remarkable polar parametriza- 
tion for the curve obtained by solving Eq. (14). Namely 
[7]: 



4a£ + f3 2 



Aan 
a 



(22) 



Once the spectrum of £ is known, Eqs. (21) and (22) 
allows one to determine the local geometry of the curve 
without need of further integration. 



B. Analytic solution in two dimensions 

In this section we will construct the analytical solution 
of the shape equation (14) and describe some its geomet- 
rical properties. The derivation presented here follows 
mostly that of Ref. [31]. As a first step we note that Eq. 
(14) has a trivial solution consisting of a circular disk 
whose radius R satisfies the cubic equation: 



aR 3 + f3R 2 - a = , 



(23) 



for all parameter values, with f3 enforcing the condition 
L = 2ttR. For large enough surface tension, however, 
we might expect this configuration to become unstable 
and the disk to buckle into a more complex shape. To 
test this hypothesis we analyze the stability of the disk 



with respect to a small periodic displacement in the radial 
direction SR = esinn^, where <p is the polar angle, n and 
integer and e a small amplitude. Expanding the energy 
(10) to second order in e yields: 



E ~ E ° + ^3 l aR3 + P n2R2 + a ( 2 ™ 4 



■ 5n 2 + 2)] , (24) 



where Eq = cf{itR 2 ) + 2tt R(a/R 2 + f3>) is the energy as- 
sociated with the circular shape. Here the Langrange 
multiplier j3 is given by Eq. (23) with R = L/2tt, so that 
j3 = [a-a(L/27r) 3 ]/(L/27r) 2 . Replacing this in Eq. (24), 
it is easy to very that the coefficient of the second order 
term becomes negative when: 



— > 167T 3 (n 2 - 1) . 
a 



(25) 



Thus the first mode that goes unstable is that associated 
with an elliptical deformation with n = 2, corresponding 
to the critical value aL 3 /a > 487T 3 . 

In order to calculate the curvature after buckling, let 
us first make Eq. (14) dimensionless by taking t = s/L, 
n = Ln, J3 = [3L 2 /a and a = crL 3 /a. We then have: 



2k" + n 3 - f3n + a = . 



(26) 



Without loss of generality, we can choose t = as the 
point where the derivative of n vanishes: n'(0) = 0. In 
addition, we can take: 



n(0) = n(l) = £ , 



(27) 



and determine the constant no later on by using the fact 
that the curve is closed, thus its integral must satisfy the 
theorem of turning tangents (4). Integrating Eq. (26) 
over n and using Eq. (27), we obtain: 

{k'f + \{k A - kt) - i/?(/c 2 - /eg) - a(k - «o) = . (28) 

Introducing the new variable y = I /(no — ft), we can 
reduce the order of the nonlinearity by one unit: 

(y') 2 = (k 3 -pko-&)y 3 + (lp-lk 2 )y 2 + k y-\. (29) 
This equation is of the form: 

(y'f=P(y)=h 2 (y-a)(y-b)(y-c), 

with h 2 = nl — j3ko — cr and a, b and c the roots of the 
cubic polynomial P(y), and is suitable to be solved in 
terms of elliptic functions. Next we make the assumption 
that P(y) has a single real root y = a and a pair of 
complex conjugate roots y = j3 ± ij. Thus P(y) admits 
the factorization: 

P(y) = h 2 (y-a)[(y-l3) 2 + 7 2 }, (30) 

which allows us to calculate the elliptic integral [32]: 






dy _i _ x (y-zi 
— = uj en , m 



(31) 



where z\ and Z2 are the roots of the quadratic equation: where II is the complete elliptic integral of the third kind: 



z 2 - 2az + 2a(3 - ((3 2 + 7 2 ) = , (32) 

and uj and m are given by: 

2 h(z 1 - z 2 ) 2 /3 - z 2 , . 

uj = , m = . (33) 

Z Z\ — Z2 

Finally, solving Eq. (31) for y and going back to our 
original variables, we have: 



n = kq 



1 — cn(co>£, m) 
z\ — Z2 cn(cjt, m) 



(34) 



In Eqs. (31) and (34) we use the standard notation for 
Jacobi elliptic function [33]. Namely, given the incom- 
plete elliptic integral of the first kind: 



u = F((j),m) 



Jo vT 



dt 



V 1 — m 2 sin 2 t 



with < m 2 < 1 the elliptic modulus, then <p is the 
Jacobi amplitude: <\> = am(w,m) and cn(^, m) = cos(/>. 
As it was noted in [31], starting from the assumption 
that polynomial P(y) has three real roots instead of one, 
would have led to a contradiction. Indeed, the fact that 
P(y) has three real roots implies that there are four dis- 
tinct values of k for which k! = 0. The solution obtained 
by assuming three real roots has however only two dis- 
tinct extrema: hence a contradiction. 

The solution (34) satisfies by construction the bound- 
ary condition k(0) = ko an d k'(0) = 0. In order for it to 
be a legitimate solution of the problem, we further need 
it to be periodic, so that k(0) = k(l) and satisfy the the- 
orem of turning tangents (4). Periodicity can be easily 
implemented by recalling that cn(x + AnK(m)) = cn(x), 
where K(m) = F(|,m) is the complete elliptic integral 
of the first kind and n an integer. This results in the 
following condition for the frequency uj: 



4nK(m) . 



(35) 



Now, the number of extrema of the function (34) in the 
interval < t < 1 is given by 2n, thus by virtue of the 
four- vertex theorem (Sec. II), n > 2. In the next section 
we will see how different values of n give shapes with a 
different number of lobes and that all of them are local 
minima of the energy (10). From the theorem of turning 
tangents (4) for a simple closed curve, on the other hand, 
we obtain: 

ko — 2tt = 



(36) 



1 

ziz 2 K(m) 


z\K(m) - 


VI -m 2 \z\ 


m 2 Y 
m 2 — 1/ 



II(n|m) = / — 
Jo (1 



dt 



(1 — n sin 2 t) v 1 — m sin 2 t 



(37) 



Solving simultaneously the transcendental equations (35) 
and (36), the quadratic equation (32) and its associated 
cubic, allows to calculate ko and j3. The solution is then 
complete. 



C. Buckled shapes 

With the help of the machinery developed in the pre- 
vious section, we can now explore the variety of shapes 
that can be accessed by the system upon chancing the 
surface tension. In Sec. IIIB we have shown that for 
any n > 2, there is a critical value <Jb n = 167r 3 (n 2 — 1) 
such that for a > <Tb n the circular shape is unstable to 
a buckled shape having n— fold rotational symmetry. A 
theorem by Tadjbakhsh and Odeh [23] guarantees that 
for any n > 2 and a > a\^ n there is a unique buckled 
solution. This is given by Eq. (34), with the parameters 
z\ and Z2 (thus uj and m) given by Eqs. (35) and (36). 
Some example of these buckled shapes are given in Fig. 
4 for various n and a. 

After buckling, increasing the surface tension results in 
a reduction of the area enclosed by the rod, with the con- 
vex regions of the curve moving away from the center of 
mass and the convex regions shifting toward the interior, 
consistently with Eq. (22). This trend continues until 
one reaches a second n— dependent critical surface ten- 
sion <j cn at which two or more points on the opposite sides 
of the boundary come into contact. cr cn is called the con- 
tact surface-tension (or contact pressure) and was first in- 
troduced by Flaherty and coworkers [24] . To gain insight 
about this regime one can consider the parametrization 
given in Eq. (21) in the center of mass frame. Symmetry 
requires that a point of contact r • n = 0, thus from Eq. 
(21) we find that the curvature at a point of contact is 
given by: 



k(t c 



/3, 



(38) 



where t c = s c /L is the rescaled arc-length coordinate of a 
point of contact and the choice of the minus comes from 
the fact that contact occurs in the concave regions of the 
rod. In the physically relevant case of n = 2, contact 
occurs at a minimum of the curvature thus: F(t c ) = 
and k(t c ) = k^. Then, using Eq. (38) and solving Eqs. 
(35) and (36) for ko and a we get ko = 2i\ and <r C 2 = 2603. 
It is worth to stress that at the point of contact, the 
curvature of the two-lobbed shape matches that of the 
circular ring before buckling. 

The same mechanism occurs for n > 3, in this case, 
however, contact does not occur at the point of minimal 
curvature, thus k(t c ) ^ ko. In order to find the contact 
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FIG. 4: Gallery of shapes obtained by solving Eq. (14) for various surface tension values and with rotational symmetry index 
2 < n < 5. The curves in the third column show the onset of point-contact at a = a cn (compare with Table I). The curves in 
the fourth column, on the other hand, show the onset of line-contact at a — ao n . In this case the curvature vanishes wherever 
the curve is in contact with itself. 



t = 



n Z 6 4 5 b 

a cn 2603.0 10740.6 25719.8 48535.5 79910.4 
& 0n 5127.7 40584.9 102795.4 203145.1 349245.2 

TABLE I: The surface tension a cn and ao n at the onset of 
point and line contact for n < 6. The units used here differ 
from those used in Refs. [8, 24]. These and can be recovered 
by dividing the values reported in the table by 167T 3 . 



surface tension a c one needs to look for a solution of the 
Frenet-Serret equations associated with the Eq. (34) that 
contains two points points of arc-length s c < s' c < L/n at 
which r(s c ) — r(s f c ) and t(s c ) = —t(s f c ). This task was 
accomplished in Ref. [23, 24] using a shooting method. 
An alternative and simpler method, was introduced by 
Djondjorov et al. [8] and consists of solving Eqs. (35) 
and (36) together with a further equation for the turn- 
ing angle at the point of contact. Some values of <r cn 
calculated with this method are tabulated in Table I. 

As it was explained by Flaherty et al. [24] and ex- 
perimentally confirmed in a beautiful work by Mora et 
al. [19], as a is increases above <r cn , the opposite sides 
remain in contact at one pair of points in each lobe. Si- 
multaneously, the stress resultant F(t c ) increases and the 
curvature decreases until it becomes zero at certain sur- 
face tension &o n > a cn . From Eq. (38) this also implies 
that P — at a — &o n . Fig. 3 shows the experimental 
realization by Mora et al. [19] of the self-contacting lobe 
together with the theoretical prediction (dashed line). 

For the simplest case of n = 2, both the critical sur- 
face tension <Jo2 and the shape of the lobes can be found 
straightforwardly by setting t — at the tip of a lobe 
(where the curvature is maximal) and t c = 1/4 (Fig. 5). 
Next we can look for a function of the form given in Eq. 
(34) that satisfies: 



A0 = 



dtk 



TT 

2 



(39) 



The second condition reflects the fact that the tangent 
vector rotates by tt in one counterclockwise loop of the 
lobe, thus it rotates by 7r/2 if we start from the mid-point 
at t = 0. Solving numerically Eqs. (39) for ko and one 
finds &o — 28.8164 and the value of &Q2 reported in Table 
I. 

For n > 3 contact does not occur at the center of the 
ring, thus finding &o n and the shape of the rod requires 
to introduce further equations to determine the position 
t c of the contact point for each lobe. As it was suggested 
in Ref. [24] this problem can be addressed by dividing 
the length of a half- lobe in two intervals: from the tip of 
the lobe to the point of contact, < t < t c , and form 
the point of contact to the end of the lobe where the 
curvature is minimal and the curve is closer to its center 
of mass, t c < t < l/n (Fig. 5). In each of these intervals 




FIG. 5: Geometry at the onset of contact. In the case of 
the 2— fold deformation mode, contact occurs in the center 
of mass of the curve, thus t c = s c /L — 1/4 if we set t — 
at the tip of a lobe, where the curvature is maximal. For 
the generic n— fold mode (with n > 3), contact occurs before 
the curvature minimum at t = l/n. Thus calculating the 
curvature requires to separate a half- lobe in two intervals, 
< t < t c and t c < t < l/n, and solve two separate problems 
with the matching condition given in Eqs. (40). 



we can then search for a solution that satisfies: 



k(t c ) = &'(0) = 
k(t c ) = k\l/n) = , 



dtk 



TT 







dtk 



TT TT 

n ~ 2 



(40a) 
(40b) 



Once again, we can use the general solution (34) to calcu- 
late the integrals, but in the second interval this should 
be transformed in such a way the derivative of the cur- 
vature vanishes at the end of the lobe, thus: t — >• 1/n — t. 
(this corresponds to translate the origin to the end of the 
lobe and invert the orientation) . Thanks to the availabil- 
ity of an analytical expression for the curvature, solving 
Eqs. (40) numerically is not hard and one can readily find 
the surface tension values reported in Table I as well as 
the shapes shown in Fig. 4 (fourth column). For a reason 
that will be explained below, all the lobes constructed in 
this way are similar regardless of n, while it is only the 
shape of the central loop that changes becoming more 
and more convex as n increases. 

Upon increasing the surface tension above ao n , the 
points of contact becomes lines of increasing length and 
the process continues until, in the limit of infinite <j, the 
shape is completely collapsed into a star consisting of n 
lines segments meeting in the center and regularly sepa- 
rated by an angle of 2tt /n. As it was observed in Refs. 
[8, 24], the shape of the partially collapsed loops can be 
constructed by rescaling that at ao n . Eq. (14) is indeed 
invariant under the following scaling transformation: 



(s, ft, /3,a) —> ( Xs 



K (3 (J 

A' A2' A3 



(41) 



Consequently, the shape of the end-loop obtained for 
a > <Ton is similar to that at &o n up to a scaling fac- 
tor A = (fron/a) 1 / 3 < 1. Accordingly, the length and 
the area of the rescaled shape are rispectivelly XL and 
A 2 A, with A the area of the reference shape obtained at 
&o n . This beautiful geometric property can be translated 
into the following algorithm to construct shapes with ex- 
tended regions of contact [8, 19, 24]: 1) given &o n and a 
surface tension a > <Jo n we construct a reference shape 
associated with &o n and calculate the scaling factor A; 2) 
we then rescale the reference curve so that its length is 
AL; 3) finally we replace each isolated point of contact 
by a line segment of length: 



£■■ 



L 

2n 



(1-A), 



(42) 



so that the total length is once again equal to L. Some 
example of shapes constructed with this procedure are 
shown in Fig. 4 (fifth column). It should be stressed 
that in the construction given here, the length of the 
contact lines is shorter by a factor 1/2 with respect to 
that given in Ref. [8]. The reason of this discrepancy 
is the following: if we think of the boundary curve as 
a physical rod, then the segments of line connecting the 
lobes with the central loop consists of two portions of 
rod in contact with each other and their length should 
be double counted. Hence the factor 1/2 included in Eq. 
(42). 

Fig. 6 shows a plot of the energy (10) rescaled by a fac- 
tor a/L for two families of solutions of the shape equation 
(14), together with some representative configurations. 
As one can see, solutions of larger n are always energeti- 
cally more expensive than the simple two-lobbed shape, 
thus are never found in a typical experimental realiza- 
tion of the problem [19]. It appears feasible, however, to 
take advantage of the initial configuration of the ring to 
bias the system toward a shape having n > 3 and pos- 
sibly use the various buckling instabilities discussed here 
to achieve conformational changes, in the same spirit of 
Ref. [34]. 




FIG. 6: The rescaled energy (10) as a function of a for the 
2— and 3— fold symmetric solution of the shape equation (14). 
The curves in the inset correspond to the values of a marked 
by the dots. The dashed- line shows the energy of the circular 
solution E /(a/L) = &/(4tt) + 47T 2 . 



structure described by the energy (10) provides a tremen- 
dous example of polymorphism and multi-stability. By 
simply increasing a single physical parameter, the surface 
tension, one can access an increasingly large number of 
stable configurations and in the limit a — »> oo the sys- 
tem admits a countably infinite number of minima cor- 
responding to fully collapsed stars having an arbitrary 
integer number of rays departing from the center. In 
addition, the problem provides a fascinating example of 
contact and cusp geometry. Finally, Eq. (14) can be 
solved in a completely analytical manner (proviso solv- 
ing numerically some transcendental equation involving 
elliptic integrals) and the solution is in excellent agree- 
ment with the data (see Fig. 3). In Sec. V we will see 
how the techniques and the intuition built up here will 
turn into an useful set of theoretical tools to address an 
important problem in cell mechanics. 



IV. THE EULER-PLATEAU PROBLEM IN 
THREE DIMENSIONS 



D. Final remarks 

If the rod is allowed to intersect, increasing the sur- 
face tension above & cn produces various self-intersecting 
curves one can also describe using Eqs. (34), (35) and 
(36). This solutions, however, do not have a physical 
relevance and for this reason they will not be discussed 
here. Interested readers will find a thorough discussion 
of this regime in Ref. [25] . 

Regardless of the applied interest to situations involv- 
ing vessels subject to a uniform pressure (thus including 
the circulatory system) , the problem reviewed in this sec- 
tion has some remarkable features. Although the physics 
involved is very simple and ultimately reduces to the 
competition between two antagonist forces that try to 
bend the rod in opposite directions, the simple composite 



In Sec. Ill we considered the problem of a soap film 
bounded by an elastic rod and constrained to lie on the 
plane. It is then very natural, and in fact tempting from a 
practical point of view, to lift the constraint of planarity 
and ask what would be the shape of a soap film bounded 
by an inextensible but flexible piece of wire that is free 
to move in three dimensions. This problem links two of 
the oldest problems in geometry and physics. Namely 
the problem of finding the shape of an Elastica subject 
to some external force and the Plateau problem of mini- 
mal surfaces. In geometry, the classical Plateau problem 
consists of finding the surface of least area that spans 
a given rigid boundary curve. A physical realization of 
the problem is obtained by dipping a stiff wire frame 
of some given shape in soapy water and then removing 
it. The original variational problem was formulated La- 
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FIG. 7: A physical realization of the Euler-Plateau problem. 
(A)-(D) Examples of minimal surfaces bounded by an elastic 
line obtained by dipping (and removing) a loop of fishing line 
(0.3 mm diameter) in soapy water. The interplay between the 
forces due to surface tension of the film and the elastic forces 
arising at the boundary can be quantified via the dimension- 
less number a — aL 3 /a with a the surface tension, L the 
length and a the bending rigidity of the boundary. For small 
values of <j, the soap film takes the form of (A) a planar disk, 
but upon increasing a it buckles in the form of (B) a saddle, 
(C) a twisted figure eight and (D) a two headed racket- like 
structure. Reproduced from Ref. [9]. 



grange [35] in the eighteenth century, however, the sub- 
ject was lifted to a whole new level by Plateau in the 
mid-nineteenth century using a series of beautiful exper- 
iments which showed that a physical realization of these 
objects naturally arises through a consideration of soap 
films [36]. Since then, the subject has inspired mathe- 
maticians [37-39] , scientists [40-42] , engineers and artists 
alike [63]. 

Since experiments with soap films are often carried 
out by dipping closed wire frames into soap and then 
pulling them out, there is a natural generalization of this 
boundary-value problem that is suggested by the follow- 
ing question: what if the soap film is not bounded by a 
rigid contour, but has instead a soft boundary such as a 
flexible, inextensible wire? The resulting problem, that 
was named Euler-Plateau problem in Ref. [9], provides 
a fascinating and highly nontrivial generalization of the 
planar problem reviewed in Sec. III. 

Fig. 7A-D, shows a gallery of shapes obtained via a 
simple experiment of dipping a loop of fishing-line into 
a solution of water and dish soap, and pulling it out. 
The resulting shapes are quite varied: short loops span 
a planar disk (Fig. 7A), intermediate loop lengths cause 
the film to spontaneously twist out of the plane (7B) and 
eventually form a planar figure-eight shape (Fig. 7C). 
At the center of the eight (where the boundary crosses) 
the surface normal rotates by 180°, and the spanning 
minimal surface is helicoid-like. Increasing the length 
even further leads to a a two-headed racket-like structure 
(Fig. 7D). 



A. Shape equation 

As in Sec. Ill we start from the simple energy given 
in Eq. (10), where M is now a mapping of a two- 
dimensional disc in 1R 3 and the boundary dM is a space 
curve. For simplicity, we assume that the region of con- 
tact between the film and the bounding elastic rod is free 
to slide, so that the film can apply normal forces, but no 
torques. From this assumptions it follows that, if the fila- 
ment is initially untwisted, it will remain untwisted even 
when buckled and bent out of the plane, although it can 
develop geometric torsion. To derive the Euler-Lagrange 
equations for the functional (10) we use again the Frenet 
frame introduced in Sec. II and consider a small virtual 
displacement of the boundary: r — > r + t;n + 776, where n 
is a displacement along the binormal direction and van- 
ishes identically if the film is constrained to lie on the 
plane. Calculating the linear variation yields [9]: 



ds (an 2 + /3) 



dM 



<b ds[a(2n" ' + k 3 -2t 2 k)-Pk\£ 

JdM 

fai ds (4k/t + 2kt')t} . (43) 



The linear variation of the area bounded by the curve 
can be expressed in a standard coordinate basis gi = 
d{R (with R position vector of the film and d{ partial 
derivative with respect to the i-th coordinate) as [43] 

S f dA= I dA(V-u- 2Hw) 
Jm Jm 

= -<£ dsT-u-2 f dAHw. (44) 

JdM Jm 

Here u and w are respectively the displacement along 
the tangent plane of the film (expressed in the basis gi) 
and its normal direction TV, H is the mean curvature of 
the film, T = N x t is the tangent-normal vector of the 
Darboux frame defined in Sec. II and the last equality 
comes from the divergence theorem on a surface. Com- 
patibility demands that the bulk and boundary variations 
must be consistent with each other. Thus, using Eq. (7) 
to express T in the Frenet frame, after simple algebraic 
manipulations we obtain: 



T • u = r\ cos $ b — £ sin # n 



(45) 



where d is the angle between the surface normal to the 
film and the normal to the boundary curve and will be 
herein referred to as contact angle. Enforcing the condi- 
tion that the energy variations (43) and (44) vanish while 
satisfying the compatibility condition (45) yields: 



1 



2a 



2a 



sin d = . 



2k' r 



H = 0. 



a 
2a 



cos $ = , 



(46a) 



(46b) 



(46c) 
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Alternatively, one can write Eqs. (46a, b) in a form that 
explicitly contains the geodesic and normal curvature of 
the boundary by expressing cos^ = n n /n and sin^ = 

KgJK. 

Unlike the planar case, the Euler-Lagrange equations 
(46) are very difficult and perhaps impossible to solve an- 
alytically for the general case. It is possible nonetheless 
to infer a number of properties without knowledge of the 
exact solution. Multiplying Eq. (46b) by k allows us to 
recast the resulting expression as (ft 2 r)' = — (a/2a)^ n , 
which integrated along the length of the closed bound- 
ary yields the following integral formula for the normal 
curvature K n : 



f • 

JdM 



ds k,„ = , 



(47) 



valid for all values of the physical parameters. Further- 
more, it is possible to consider a special class of solutions 
of Eqs. (46) for which the contact angle i? is constant 
along the curve. Then Eq. (47) can be used to prove that 
$ must then be 90°. To understand this latter statement 
we first need to recall that the normal vector n, thus 
the contact angle #, is undefined at inflection points (i.e. 
where n = 0, see Sec. II). In particular, if d is constant, n 
must be everywhere positive. Then, writing n n = n cos d 
in Eq. (47) yields: 



cost? 



ds k = , 



dM 



(48) 



but since n is strictly positive, this last condition is pos- 
sible if and only if cos$ = 0, thus d = ±tt/2: i.e. K n = 
and the surface normal must be perpendicular to the 
curve normal everywhere if d is to be constant. This 
imposes strong constraints on the nature of the solution 
in this case as we will now argue. 

Curves having zero normal curvature everywhere 
are called asymptotic curves and their tangent vectors 
asymptotic directions. Every curve that lies on the plane 
is clearly asymptotic. On a minimal surface, on the other 
hand, there are two asymptotic curves passing through 
each point where the Gaussian curvature K is non-zero 
and these curves always intersect at 90°. If the min- 
imal surface contains flat points (i.e. isolated points 
where K — 0), these are traversed by JV > 2 asymptotic 
curves crossing at n/Af and forming, in their neighbor- 
hood, n valleys separated by ridges [44]. These obser- 
vations suggest that the only simply connected minimal 
disk bounded by a closed asymptotic curve lies in the 
plane; otherwise the asymptotic directions in the disk 
would form a vector field tangent to the boundary. In 
a simply connected disk this vector must then enclose 
a singularity or vortex, in the neighborhood of which 
the structure of the asymptotic directions would disagree 
with the previous classification. Hence, we conclude that 
the only solution having constant contact angle is given 
by the flat disk. 



B. Numerical simulations of the soap film shapes 

To explore the variety of possible shapes resulting from 
the solution of Eqs. (46) and the transitions between 
them, Mahadevan and the author [9] minimized a dis- 
crete analog of the total energy (10). The soap film is 
approximated as simplicial complex consisting of an un- 
structured triangular mesh. The internal edges of the 
triangles are treated as elastic springs of zero rest-length 
so that the total energy of the mesh is given by: 



E = a Y, (s v ) K 2 v + kY^\e? 



E 

vedM 



(49) 



eeM 



The first sum runs over the boundary vertices and (s v ) = 
(s v + s v -i)/2 is the average of the length of the two edges 
meeting at v. The curvature of the boundary is calcu- 
lated as k v = \t v — t v -i\/(s v ) with t v and t v -\ the tan- 
gent vectors at v. The second sum in Eq. (49) runs over 
all internal edges. If the triangles are equilateral, this 
yields a discrete approximation for the soap-film energy 
with the spring stiffness proportional to the surface ten- 
sion, i.e. a ~ 4/c\/3/(2 — B/E), where B/E is the ratio 
between the number of boundary edges B and the total 
number of edges E of the triangular mesh. The choice 
of minimizing the squared length of the edges, instead 
of the area of the triangles, is motived by numerical sta- 
bility. Replacing the spring energy in (49) with the sum 
of the area of the triangles, has the effect of shortening 
the range of interaction between the boundary and the 
interior to the single strip of triangles at the boundary. 
Most standard local optimization algorithms would then 
attempt to reduce the area of these triangles to zero, 
thus suppressing the interaction between the boundary 
and the interior. In the continuum limit, the energy (49) 
approaches that given by (10) for original problem, and 
the sequence of shapes obtained with this method is in 
excellent agreement with our experimental observations 
shown in Fig. 7. 

Numerical minimization of Eq. (49) using a conju- 
gate gradient method leads to a variety of shapes de- 
pending on the values of a, /c, the boundary length L 
as well as the initial shape of the domain. An interac- 
tive gallery of the shapes obtained from the numerical 
simulations is available on-line [64]. A first set of sim- 
ulations were run by starting with an unstructured tri- 
angular mesh of 474 vertices bounded by an elongated 
hexagonal polygon; the mesh points are started out ran- 
domly displaced transverse to the plane. As the spring 
constant k was increased, corresponding to an increase 
in the surface tension, the system moved through a se- 
ries of transitions shown in Fig. 8. Beyond a critical 
value of the surface tension (relative to the bending stiff- 
ness, as discussed earlier), the planar disk buckled into a 
two-fold mode consistent with the linear stability analysis 
reported in Sec. IIIB. For kL 3 /a « 643, the planar ellip- 
tical shape transitions to a twisted non-planar saddle-like 
shape, as experimentally observed. This is accompanied 
by an increase in the normal curvature of the boundary 
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and there is a progressive twisting of the central "waist" 
of the surface. A suitable quantity to characterize the 
non-planarity of the boundary is obtained by integrating 
its absolute normal curvature: 



(\K n \) = (b ds\K n \ 

IdM 



(50) 



A plot of (|tt n |) for a wide range of kL 3 /a values is shown 
in Fig. 8. A further increase in k leads to an increase in 
the twist until the surface normal undergoes a full 180° 
rotation across the waist for kL 3 /a « 740. The surface 
then reverts to a planar configuration in the shape of the 
number eight (see Fig. 8E). 

Once the surface reaches the planar figure-eight con- 
formation an increase in the surface tension does not pro- 
duce further conformational changes until kL 3 /a « 2000. 
Beyond this value, the surface again becomes non-planar, 
self-intersects at the pinch of the eight, and the two lat- 
eral lobes start bending toward each other. The curva- 
ture of the lobes increases with increasing k until they 
intersect to produce the complex shape shown Fig. 8H. 
The transition to this latter surface is analogous to the 
transition from the figure-eight to the two headed racket- 
like structure observed in the experimental realization of 
the problem. However, in the experimental system, self- 
intersection cannot occur and adhesion favors the forma- 
tion of a line of contact between the two lobes (see Fig. 
7D). 

Not surprisingly, given the nonlinear nature of the gov- 
erning equations (46), the solutions are not expected to 
be unique. Numerically, one finds that the final equilib- 
rium state strongly depends on the initial configuration 
of the system. If the initial configuration is not elongated, 
as in say a triangular mesh bounded by a hexagon, which 
is closer to being circularly symmetric, other possible sta- 
ble configurations are found. As in the previous case, for 
small values of the surface tension, the system rapidly re- 
laxes into a flat circular disk. However, for kL 3 /a ~ 855 
the system transitions to a saddle-like configuration that 
is the classical Enneper minimal surface. Upon increas- 
ing the surface tension still further, the curvature of the 
boundary and film becomes larger still and the surfaces 
self-intersect leading to the beautiful structure shown in 
Fig. 9D. 



C. Asymptotic analysis 

The numerical data show that the transition from 
the planar two-fold symmetric shape to the non-planar 
twisted eight is consistent with a supercritical pitchfork 
bifurcation, with (|ft n |) ~ (a — at) 1 / 2 , where the angular 
brackets stand for an average along the boundary curve. 
Indeed Fig. 10 shows a best fit of the numerical data at 
the onset of the transition, the exponent obtained from 
the fit is 0.489 ± 0.005. 

To understand the origin of the pitchfork bifurcation, 
it is helpful to introduce a specific approximate represen- 



tation of the twisted soap film in terms of the following 
one-parameter family of surfaces: 



r(l 



t 2 ) COS( 



y = r (1 — t 2 ) sin0, 
z = £(r 2 /#)sin20, 



(51a) 
(51b) 
(51c) 



where < < R and < < 2tt are the usual plane- 
polar coordinates. Here — 1 < t < 1 is a parameter 
that characterizes the family of shapes. Eq. (51a) is 
the simplest parameterization of a twisted saddle with 
disk topology and whose degree of non-planarity is con- 
trolled by the parameter t. For t — 0, the surface given 
by (51a) is a disk of radius R; as t becomes positive (neg- 
ative), the surface becomes a right-handed (left-handed) 
twisted figure-eight, and for t = ±1 the surface reduces 
to a flat twisted figure-eight bounded by a lemniscate of 
Gerono. In the following we will exploit the topological 
equivalence between this family of surfaces and the ac- 
tual twisted soap film to study the global properties of 
the surface at the onset of the transition, when \t\ <C 1. 
Lengthy, but straightforward calculations give an ex- 
pression for the curvatures of the boundary up the fourth 
order in t: 

2b 2t 3 

K<n ~ — — sin 20 -\ — — (3 sin 20 — sin 40 + sin 60) , 



R ^ R 

1 t 2 



R R 



(1 + 3 cos 20- 5 cos 40) 



Similarly, the arc-length is given by: 

ds 2 = # 2 [(l + £ 2 ) 2 -2£ 2 (cos20-cos40)] . (52) 

From this one can calculate an approximated expression 
for the total energy at the onset of twisting (i.e. for 
t~0). Namely: 
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where Eq is again the energy of the circular configuration. 
Finally, the condition of inextensibility for the boundary 
curve is given by: 

L = <L dsxj 2ttR [1 + t 2 (1 - \t 2 )] . (53) 

JdM 

Using equation Eq. (53) to eliminate R and then taking 
the derivative of the approximate energy with respect 
to t yields the following equation of equilibrium for the 
amplitude parameter t that characterizes the shape of 
the boundary: 



£(96tt 3 - a) + | at 3 = , 



(54) 



which is the normal form of a supercritical pitchfork bi- 
furcation and implies: 



(K|) ~ (a - at) 1 / 2 , a t = 967r 3 



(55) 
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FIG. 8: Some solutions of the Euler-Plateau problem and a bifurcation diagram. Shapes obtained by minimizing the discrete 
approximation (49) of the continuum energy of a minimal surface bounded by an elastic line (10) as a function of the dimen- 
sionless number kL 3 /a. For small values of the parameter, the flat circular disk of radius L/2tt has the lowest energy. Upon 
increasing kL 3 /a the disk first buckles into flat elliptical shape (see text). For kL 3 /a > 643, this elliptical disk deforms into a 
three-dimensional saddle-like shape, becoming increasingly twisted with an increase in kL 3 /a > 643 until the surface eventually 
becomes flat again, but now forming a fully twisted eight-shaped conformation. This twisted flat shape is a local energetic 
minimum for 740 < kL 3 /a < 2000. For still larger kL 3 /a the twisted surface self-intersects at the waist of the figure-eight and 
the two lateral lobes bend toward each other, while for very large kL 3 /a the surface exhibits several self-intersections leading 
to the complex structure (H). Below the shapes we show the bifurcation diagram for the system, characterizing the shape using 
the absolute normal curvature of the boundary integrated along its length (a measure of the amplitude of the instability) as 
a function of the bifurcation parameter kL 3 /a, along with the location of the transitions described above. Reproduced from 
Ref. [9]. 



having used that (|ft n |) ~ t. Analogously one can cal- 
culate the behavior of the Gaussian curvature of the 
bounded minimal surface at the onset of twisting. Us- 
ing the Gauss-Bonnet [14] one has: 



from where, the integrated Gaussian curvature (K) 
J dA K is given by: 



(K) ~ (a - & t ) 



(57) 



/• 

J M 



dAK = 2tt 



ds K n « Ant 



dM 



(56) 



in excellent agreement with the numerical simulations 
shown in Fig. 10. 
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FIG. 9: Further equilibrium solutions of the Euler-Plateau problem. An initial configuration of a triangular mesh bounded by a 
hexagonal perimeter gives the classic Enneper minimal surface rather than the figure-eight shape shown in Fig. 8. Reproduced 
fromRef. [9]. 



D. Final remarks 

The Euler-Plateau problem, reviewed in this section, 
is a fascinating example of a free boundary value prob- 
lem obtained from the natural marriage of two classic 
problems of geometry and physics. While a minimal 
physical realization of the Euler-Plateau problem is a 
kitchen-sink experiment, its ramifications are likely to 
go far beyond this specific manifestation, just as the me- 
chanics of soap films and elastic filaments have been rel- 
evant for the study of matter, not just at the every-day 
scale, but also for systems that range from molecules [40] 
to black holes [42]. However, unlike its planar version 
reviewed in Sec. Ill, the three-dimensional problem ap- 
pears as a formidable mathematical challenge and the 
open questions outnumber the answers by far. Some im- 
mediate mathematical questions include proofs of exis- 
tence, regularity of solutions, as well as measures of the 
non-uniqueness of solutions, particularly for large values 
of the only parameter in the problem a or its discrete 
analogue kL 3 /a. 



V. APPLICATION TO CELL MECHANICS 

As we pointed out in the previous sections, determin- 
ing the optimal shape of a soap film bounded by an elas- 
tic element is not just an intriguing mechanical puzzle 
inspired by the observation of the mundane, but an im- 
portant prototype problem for various physical systems 
whose geometry is softly constrained. In this section we 
will consider an important application of these methods 
to the case of contractile cells adhering to compliant sub- 
strates. 

The propagation of mechanical forces generated dur- 
ing cell-matrix adhesion is now recognized as one of the 
fundamental mechanism beyond a variety of cellular pro- 
cesses. By mean of adhesive contacts with the substrate, 
a cell is able to sense mechanical cues from the environ- 
ment, transmit them to its interior and transduce them 



into chemical signals that in turn activate a spectrum 
of biological responses [45-47] . For instance cells adher- 
ing to softer substrates spread less and prefer to have 
well rounded morphologies, while they are more likely 
to exhibit branched patterns on stiffer substrates with 
greater spread area [48, 49] (Fig. 11). Experiments on 
micro-patterned adhesive islands revealed that cell pro- 
liferation and spreading sensitively depend on adhesion 
geometry [50] . Perhaps the most impressive effect of the 
cell-substrate interaction is its ability to determine the 
fate of some type of stem cells that, adapting themselves 
to the substrate stiffness, might becomes either neurons 
or bone cells [51]. 



A. The contractile film model for adherent cells 

Fully spread cells adhering on a solid substrate gen- 
erally exhibit a flat morphology with the nucleus form- 
ing bump in the center. Adhesion with the substrate 
predominately occurs along the cell boundary in a dis- 
crete set of adhesion sites or continuously. This suggests 
that insight about the geometry and the mechanics of 
cell adhesion can be gained, in first approximation, by 
treating the cell as a two-dimensional film that is in con- 
tact with the substrate through points localized along its 
contour. At time scales when the cell is fully spread, 
the forces generated by the actomyosin cytoskeleton are 
mainly contractile and this gives rise to an effective ten- 
sion that is transmitted to the substrate thorough the 
focal adhesions [52]. The actin cortex localized at the 
cell periphery naturally resists to the inward contraction 
by opposing an elastic force. The overall effect of ac- 
tomyosin contractily, cortex elasticity and adhesion on 
the shape of the cell was investigated in Ref. [13], by 
mean of a simple and yet very rich mechanical model in- 
spired by the planar film problem reviewed in Sec. Ill: 
the Contractile Film Model. In this model, the shape of 
the cell-substrate contact line is parametrized by closed 
plane curve and the total mechanical energy of the system 



15 




647 648 649 650 651 652 653 

^3 , 



kL 6 /a 



-0.05 



(K) 



-0.15 




■0.10 



kL 6 /a 



FIG. 10: Characterizing the onset of twist and non-planarity. 
The left panel shows the average absolute normal curvature 
of the boundary (|ft n |) as a function of kL 3 /a in the neigh- 
borhood of the onset of the transition to non-planarity. The 
solid line on the left figure is obtained from a best fit of the 
data. The exponent resulting from the fit is 0.489 ± 0.005, 
and is consistent with the instability being of a supercritical 
pitchfork type. The right panel shows the average Gaussian 
curvature of the surface (K) as a function of kL 3 /a in the 
neighborhood of the onset of the transition to non-planarity. 
Reproduced from Ref. [9]. 



is approximated, on the basis of symmetry arguments, in 
the following form: 

E = a dA + (b ds(aK 2 + (3 + k s p\r-r \ 2 ). (58) 

JM JdM 

Here a represents the effective surface tension in the cell 
due to cytoskeletal contractility, while the last term in 
Eq. (58) represents the strain energy induced by the cell 
on a substrate of stiffness k s through focal adhesions lo- 
calized at the cell edge [53] with density p(s), so that the 
total number of adhesions is Na = § dsp. The reference 
configuration ro(s) describes the shape attained by the 
cell boundary on a rigid substrate, and is a constraint set 
by the adhesion geometry. It is important to stress that 
all the physical constants in Eq. (58) are effective as, in 
addition to the normal mechanical properties of the pas- 
sive elements (the cell membrane, the liquid component 
of the cytosol etc.), they account for the additional forces 
due to the actomyosin activity. 

A class of theoretical models by Bischofs et al. [11, 12] 
have been constructed around the competition of bulk 




FIG. 11: Fluorescently labeled actin cytoskeleton in fibrob- 
lasts adhering on a soft (left) and stiff (right) substrate. On 
the soft substrate the cell has a modest spread area and a 
rounded morphology, while on the stiffer substrate the cell 
spreads more and its geometry is concave and branched. Re- 
produced with permission from Ref. [45]. 



and peripheral contractility and ignored the bending elas- 
ticity of the actin cortex (i.e. a = 0). In analogy with 
the Laplace law of capillarity, the steady state cell con- 
tour is then described by concave circular arcs of radius 
f3/(j connecting adhesion sites. The model introduced in 
Ref. [13] and reviewed here, instead, focuses on the op- 
posite limit and considers the regime in which the force 
balance is dominated by the competition between cortex 
elasticity and bulk contractility, while the effect of pe- 
ripheral contractility is negligible (i.e. /3 = 0). In this 
scenario, the curvature is generally no n- uniform, espe- 
cially in the neighborhood of adhesion sites. As we will 
see in Sec. VC, incorporating bending elasticity leads to 
an extremely rich polymorphism and allows for a tran- 
sition from purely convex to purely concave cell shape 
reminiscent of that observed in experiments on cardiac 
myocytes [49]. From a purely mechanical point of view, 
the Contractile Film Model defined by Eq. (58) is equiv- 
alent to a fluid film bounded by a planar Elastica and 
attached to a continuous or discrete set of links to a com- 
pliant substrate. Unlike the problem reviewed in Sec. Ill, 
however, Eq. (58) does not involve any hard constraint 
on the perimeter of the cell, which is then only softly 
constraint by the adhesion with the substrate. 



B. Continuous adhesion 

If the periphery of the cell forms contact with a single 
continuous adhesion site, so that p — 1/L with L the 
perimeter of the cell, the energy minimization problem 
admits a simple solution. In the presence of a uniform 
and isotropic substrate, we can assume the reference con- 
figuration to be a circle of radius Ro so that a natural 
minimizer of the energy (58) would be a circle or radius 
R. Thus, setting /3 = 0, k = R~ x and p~ x = 2ttR and 
minimizing Eq. (58) yields the following cubic equation: 
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FIG. 12: Relative cell size R/Ro as a function of substrate 
stiffness k s (solid black circles) for smooth muscle cells, 4 
hours after plating on continuous elastic gels [54] . Cell radius 
is estimated from the projected cell area reported in Ref. [54] 
as R = v / area/7r. Substrate stiffness k s is determined from 
substrate Young's modulus E s as : k s = aE s , where a is the 
characteristic focal adhesion size, with a ~ 1 \±m. Solid (red) 
line represents the solution to Eq. (59) with with a = 1.05 
nN/jLim and a/Ro — 0.16 nN//xm. Reproduced form Ref. 
[13]. 



(k s + ira)R 3 - k a R R 2 - ira = . 



(59) 



The equation contains two length scales, Ro and £ = 
(a/a) 1 / 3 , and a dimensionless control parameter k s /a 
expressing the relative amount of adhesion and contrac- 
tion. For very soft anchoring k s <C a and Eq. (59) admits 
the solution R = £, whereas if the cell is rigidly pinned 
at adhesion sites, k s ^> a and R ^ Ro. For interme- 
diate values of k s /a the optimal radius i? interpolates 
between £ and i?o and is an increasing function of the 
substrate stiffness k s , in case £ < it^? or a decreasing 
function if £ > i?o- For £ = Rq, the lower and upper 
bound coincide, and the solution is R = Ro. In par- 
ticular, the case Ro > £ reproduces the experimentally 
observed trend that cell projected area increases with in- 
creasing substrate stiffness before reaching a plateau at 
higher stiffnesses [48, 49]. 

Fig. 12 shows a fit to the measured projected areas of 
smooth muscle cells (SMCs) adhering to continuous elas- 
tic gels of varying substrate elastic modulus [54] with the 
solution of Eq. 59. The data for the spread area of SMCs 
are taken 4 hours after plating onto the substrate, when 
they retain rounded morphologies. The fitted value for 
surface tension a = 1.05 nN//im comes to the same or- 
der of magnitude as reported for endothelial cells [12, 55], 
epithelial cells [56] and is consistent with the numerical 
estimate provided earlier. The fit also provides a value 
for the bending rigidity a = 4.62 x 10 -16 Nm 2 . The 
asymptotic behavior and various limits of the solution 
are well captured by the interpolation formula: 
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FIG. 13: Cell anchored onto three pointwise adhesions lo- 
cated at the vertices of an equilateral triangle. The curvature 
(A) and the tangent angle (B) as function of arc-length for 
crRl/a = 10, k s Rl = 50 and N A = 3. The circles are ob- 
tained from a numerical minimization of a discrete version 
of the energy (10), while the solid lines corresponds to our 
analytical approximation. (C) The total cell length £ as a 
function of adhesion stiffness. For small stiffnesses the cell 
boundary form a curve of constant width (lower inset) and 
C = 7vw, with w the width of the curve. This property breaks 
down for larger stiffnesses when inflection points develops (up- 
per inset). (D) The curvature ko at the adhesion points as 
a function of the substrate stiffness for various contractility 
values. The points are obtained from numerical simulations 
while the solid lines correspond to our analytical approxima- 
tion. Reproduced from Ref. [13]. 



indicating that larger surface tension, hence larger cell 
contractility a leads to lesser spread area, consistent with 
the experimental observation that myosin-II activity re- 
tards the spreading of cells [57]. Standard stability anal- 
ysis of this solution under a small periodic perturbation 
in the cell radius shows that the circular shape is always 
stable for any values of the parameters a, k s and Rq. 



C. Discrete adhesions 

For cells adhering to discrete number of adhesion sites, 
one can show that the circular solution for the cell bound- 
ary is never stable and there is always a non-circular 
configuration with lower energy. For simplicity, we as- 
sume that Na adhesion sites are located at the vertices 
of a regular polygon of circumradius Ro, with density 
P( s ) = 5^ n =cT d(s — it), and £ the distance between sub- 
sequent adhesions. Following the derivation given in Sec. 
Ill A and taking into account that 5(<f ds p) = 0, since 
the number Na of adhesions is assumed constant, the 
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FIG. 14: (a) a < <7 c i, cell contour is everywhere convex with constant width, (b) a = a p , cell contour is purely concave with 
cusps at adhesion points and without protrusions, (c) a > a C 2, cusps are connected to the substrate by means of a protrusion 
of length £. 



Euler-Lagrange equation can be obtained in the form: 



adhesion sites : 



a (2k" 



n a -i 
2k s J2 5 ( s 

2 = 



(r — ro ) • n = 



( 61 ) 

Due to the TV^-fold symmetry of the adhesion sites, ad- 
hesion springs stretch by an equal amount A in the 
direction of the normal vector: (ri — r^) • n^ = A, 
i = 1,2... Na- As a consequence of the localized ad- 
hesion forces, the curvature is non- analytical at the ad- 
hesion points. Integrating Eq. (61) along an infinitesimal 
neighborhood of a generic adhesion point z, one finds: 



k s 
2a 



A. 



(62) 



The local curvature of the segment lying between ad- 
hesion points is on the other hand determined by Eq. 
(14), with j3 — and boundary conditions: n(i£) = 
k((z + 1)£) = tto- Without loss of generality we con- 
sider a segment located in < s < £. Although the exact 
analytic solution this nonlinear equation is available (see 
Sec. Ill), an excellent approximation suitable for this 
purpose can be obtained by neglecting the cubic nonlin- 
earity (Fig. 12C,D). With this simplification, Eq. (61) 
admits a simple solution of the form: 



k(s) = tto + -— s(s - £) . 
4a 



(63) 



Eqs. (63) and (62) immediately allow us to derive a con- 
dition on the cell perimeter: £ = 2k s A/a. Furthermore, 
the latter condition leads to a linear relation between 
traction force T = 2/c s A, and cell size : 



T = a£ 



(64) 



which is indeed observed in traction force measurements 
on large epithelial cells [56]. 

To determine the end-point curvature ^o? we use the 
turning tangents theorem for a simple closed curve (4). 
This leads to following relation between local curvature 
and segment length, or equivalent ly traction force, at the 



K 



24a 



2tt 



(65) 



Finally, to determine the optimal length of the cell 
segment £, we are going to make use of a remarkable geo- 
metrical property of the curve obtained from the solution 
of Eq. (61) with discrete adhesions: the fact of being a 
curve of constant width [16]. The width of a curve is the 
distance between the uppermost and lowermost points on 
the curve (see lower inset of Fig. 12B). In general, such 
a distance depends on how the curve is oriented. There 
is however a special class of curves, where the width is 
the same regardless of their orientation. The simplest 
example of a curve of constant width is clearly a cir- 
cle, in which case the width coincides with the diameter. 
A fundamental property of curves of constant width is 
given by the Barbier's theorem [16], which asserts that 
the perimeter L of any curve of constant width is equal 
to width w multiplied by tt: L = ttw. As illustrated in 
Fig. 12D, this is confirmed by numerical simulations for 
low to intermediate values for contractility and stiffness. 
With our setting, the cell width is given by: 

w = (R - A)(l + costt/7V a ) + h(£/2) , (66) 

where h(s) = J ds' sin^(s / ) is the height of the curve 
above a straight line connecting subsequent adhesions 
and 6 the turning angle. For small deflections h can be 
approximated as: 



h(s) &s(£-s) 



N A £ 48a 



s(£- 



(67) 



Using this together with Eq. (66) and the Barbier's the- 
orem with L = Na^ allow us to obtain a quartic equation 
for the cell length: 



N A £ 



(I + costt/TVa) [Rq 
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FIG. 15: Traction force as a function of substrate stiffness (A) 
and contractility (B) obtained from a numerical minimization 
of a discrete analog of Eq. (10). For convenience the various 
quantities have been rescaled as follows: L — L/Ro, T = 
TR^/a, & — aRo/a and k s = kR^/a. The solid curves denote 
the approximate traction values obtained from Eq. (68). (C) 
Boundary length L obtained by increasing (squares) and then 
decreasing (triangles) the contractility for substrate stiffnesses 
k s — 100 (green squares, black triangles) and k s = 120 (red 
squares, blue triangles). The diagram shows bistability in the 
range <7 P < a < crj. (D) The critical contractility g\ and <j\ 
as functions of substrate stiffness. 



Fig. 15A,B show plots of the traction (64) with £ deter- 
mined by solving Eq. (68). The plots support the ex- 
perimental trend that traction increases monotonically 
with substrate stiffness k s before plateauing to a finite 
value for higher stiffnesses [58, 59]. The plateau value 
increases with increasing contractility (Fig. 15A). Trac- 
tion force grows linearly with increasing contractility for 
aR^/a <C 1, before saturating to the value 2k s Ro at large 
contractility aR^/a ^> 1, as shown in Fig. 15B. 



D. Inflections, cusps and protrusions 

For low to intermediate values of a and k s , cell shape is 
convex and has constant width. Upon increasing a above 
a k s — dependent threshold cr*, however, the cell boundary 
becomes inflected (see Fig. 16 and upper inset of Fig. 
12B). Initially a region of negative curvature develops in 
proximity of the mid point between two adhesions, but 
as the surface tension is further increased, the size of this 
region grows until positive curvature is preserved only 
in a small neighborhood of the adhesion points. Due to 
the presence of local concavities, the cell boundary is no 
longer a curve of constant width. 

Upon increasing a above a further threshold value 
<J2, the inflected shape collapses giving rise to the star- 
shaped configurations shown in upper right corner of 




FIG. 16: Phase diagram in plane of a — aR^/a and k s — 
ksRo/a showing optimal configuration obtained by numerical 
minimization of the energy (10) for Na — 3. 



Fig. 16. These purely concave configurations are made 
by arcs whose ends meet in a cusp. The cusp is then 
connected to the substrate by a protrusion consisting 
of a straight segment of that extends until the adhe- 
sion point rest position, so that A « (Fig. 14C). 
The cell boundary becomes pinned at adhesion sites as 
a result of having to satisfy force-balance, Eq. (61), and 
adhesion-induced boundary condition, Eq. (62), while ac- 
commodating large contractile tensions at its neighbour- 
hood. This results in spontaneous expansion in the cell 
perimeter. Unlike the previous transition from convex to 
non-convex shapes, this second transition occurs discon- 
tinuously and is accompanied by a region of bistability in 
the range <r p < a < a\. This is clearly visible in the hys- 
teresis diagram in Fig. 15C showing the optimal length 
obtained by numerically minimizing a discrete analog of 
Eq. (10) in a cycle and using as initial configuration 
the output of the previous minimization. The onset of 
bistability is regulated by substrate stiffness as shown in 
Fig. 15D, with stiffer substrates promoting transition to 
cusps at lower a\. Away from the protrusion, the cur- 
vature has still the form given in Eq. (63), with k,o = 
so that the boundary is everywhere concave or flat and 
the bending moment M = Iolkz does not experience any 
unphysical discontinuity at the protrusions origin. 

At a = <j p the protrusion have zero length and the 
corresponding configuration of the cell boundary form a 
special curve from which all the shapes having non-zero 
protrusion length can be constructed by means of the 
similarity transformation (41) discussed in Sec. IIIC. To 
see this let £ p be the protrusion length and let us set 
£ p = at a = cr p , so that the shape of the cell will 
be of the kind illustrated in Fig. 14B. Then given the 
surface tension a > <j p we calculate the scaling factor 
A = (crp/a) 1 / 3 , we rescale the reference curve by A and fi- 
nally we fill the distance between the adhesion points and 
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the cusps with straight segments of length £ p = Rq(1 — A) 
(since Ro is the circumradius of the reference shape and 
XRo that of the rescaled shape). This latter step, ulti- 
mately allows us to formulate a scaling law for the length 
of protrusions that can be tested in experiments: 



i P /Ro = 1 - (<t p /<t) 



1/3 



(69) 



This transition from a smooth shape to a self-contacting 
shape with cusps is reminiscent of the post-buckling sce- 
nario of an elastic ring subject to a uniform pressure re- 
viewed in Sec. Ill, but unlike this case, where the system 
undergoes a continuous transition from a simple curve 
to a curve with lines of contact, here the transition is 
discontinuous along both the loading branch (increasing 
k s ) and the unloading branch (decreasing k s ). The tran- 
sition has moreover a strong topological character since 
it involves a jump in the rotational index of the curve, 
whose total curvature after the transition becomes: 

dsn = tt(2- N a ) , 

dM 

by virtue of Eq. (5). Some further detail about the ge- 
ometry of protrusions in this model can be found in Ref. 
[13]. It should be stressed, however, that our knowledge 
of this phenomenon is still very preliminary and the ori- 
gin itself of this exotic instability (which in same extent is 
reminiscent of the sulcification instability in neo-Hookean 
materials [60-62]) still unclear. One of the fundamen- 
tal aspect that distinguishes the Contractile Film Model 
form classical elasticity relies on the fact that also the 
perimeter is not hardly constrained, but only subject to 
a soft constraint by mean of the adhesion springs. The 
length of an elastic object affects its overall flexibility (i.e. 
long filaments are floppy and easy to bend, while short fil- 
aments are stiff), thus, when the effective surface tension 
is increased, the whole cell boundary becomes shorter 
and stiffer. Because stiff materials are difficult to bend, 
but easy to break, a possible interpretation could be the 
following. For sufficiently large adhesion, increasing the 
surface tension has the effect of bending and stiffening 
the cell boundary in proximity of the adhesion sites, un- 
til, above a certain surface tension, the cell boundary is 
too stiff to continue bending and fractures. The cracks 
are localized at the adhesion points, where the curvature 



initially focuses, giving rise to the cusps observed in the 
simulations. However, a thorough understanding of this 
phenomenon remains a challenge for the future. 



VI. CONCLUSIONS 

In this article we reviewed some of the fundamen- 
tal mathematics and physics of softly constrained films. 
Starting from the original 19th century query of finding 
the shape of a two-dimensional elastic ring subject to a 
uniform pressure [6] , we retraced the developments of this 
concept until the most recent analytical and experimen- 
tal progress [7-9, 19] and we examined an application to 
the important problem of cell adhesion. While inspired 
by the observation of the mundane, the subject has re- 
ceived renewed interest by a heterogeneous community 
of physicists, applied mathematicians and engineers and 
is now a powerful prototype problem with potential ap- 
plications to various physical systems whose geometry 
is softly constrained. Applications to biological tissues 
and cells are by far the most promising. Although it 
relies only on a few fundamental physical mechanisms, 
compared to the extraordinarily complex orchestration 
of the cell machinery, the Contractile Film Model has 
provided some interesting insight on the geometry of ad- 
herent cells, including the experimentally observed tran- 
sition from convex to concave geometry through the for- 
mation of cusps and protrusions. While a thorough un- 
derstanding of the latter phenomena requires a deeper in- 
vestigation, the scenario suggested by these preliminary 
results is very promising: is a whole new class of biologi- 
cally relevant instabilities and morphological transitions 
lurking beyond geometrical constraints? 
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